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Autonomous Navigation Error Propagation Assessment 
for Lunar Surface Mobility Applications 

Bryan W. Welch and Joseph W. Connolly 
National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Abstract 

The NASA Vision for Space Exploration is focused on the return of astronauts to the Moon (ref. 1). 
While navigation systems have already been proven in the Apollo missions to the moon, the current 
exploration campaign will involve more extensive and extended missions requiring new concepts for 
lunar navigation. In this document, the results of an autonomous navigation error propagation assessment 
are provided. The analysis is intended to be the baseline error propagation analysis for which Earth-based 
and Lunar-based radiometric data are added to compare these different architecture schemes, and quantify 
the benefits of an integrated approach, in how they can handle lunar surface mobility applications when 
near the Lunar South pole or on the Lunar Larside. 

1. Introduction 

In support of NASA’s Vision for Space Exploration (ref. 1), error propagation of an autonomous 
navigation system is analyzed for multiple lunar surface mobility applications. This error propagation is 
conducted using Gaussian random errors as the errors associated with accelerometers and gyroscopes. In 
the Apollo era, lunar rovers used simple autonomous navigation systems to traverse the surface and return 
to the lunar module. The analysis provided in this document looks at two lunar surface mobility 
applications that model mobility applications, and determine the maximum error propagation that 
occurred due to the Gaussian random errors in the autonomous navigation system. This analysis is to be 
used as a starting point for future analysis which will include radiometric data (either from Earth-based 
assets or a Lunar-based constellation) to improve navigation performance. 


2. Mission Design 

Two lunar surface tracks are designed for the purposes of this analysis. The two tracks are quite 
different from each other in two ways: velocity profile and track route. Typically, lunar surface roving can 
be characterized by several stages such as: 

1 . Acceleration from lunar module 

2. Traversing lunar surface to destination 

3. Deceleration to stop at destination 

4. Acceleration to lunar module 

5. Traversing lunar surface to lunar module 

6. Deceleration to stop at module 

The first of the two surface roving applications takes a more simplistic approach to this surface roving 
profile, in that it eliminates the acceleration and deceleration profiles and keeps a constant velocity profile 
throughout the simulation. 
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Circular Constant Velocity Track 


The first track that is created follows the outline of a small circle centered on the moon’s fixed axis 
latitude/longitude point (0° N. 0 °E.). The radius of the circular track is 15 km on the lunar surface. 
Therefore, the rover would be at a maximum of 30 km away from the lunar module when at 1 80° rotation 
from the lunar module. The duration of the roving mission was set at 6 hr. The starting point of the track 
is placed at the point (0.496 °N. 0 °E.), corresponding to a 15 km radius. 

This circular track is designed to have a constant velocity throughout the simulation. Since the track 
has a radius of 15 km and it takes 6 hr to traverse the track, the resulting speed the rover has is 4.375 m/s. 

It is important to note that this constant speed track does not start with zero magnitude velocity, nor 
does it slow down, as a normal rover application would. This rover has a constant speed throughout the 
entire duration of the simulation. The only acceleration resulting from this profile is due to the track being 
circular. 

Figure 1 provides an illustration of the circular track, which is created in Satellite Tool Kit (STK). 
Figure 1 illustrates the constant speed of the rover on this track by capturing images at 0, 0.75, 1.5, 

2.25, 3, 3.75, 4.5, 5.25, and 6 hr after the starting point of the rover from the lunar module. Note that in 
the image bottom right region, representing the rover 6 hr after departure from the lunar module, the rover 
is again at the starting point. 



Figure 1 . — Circular track velocity illustration. 
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Straight Variable Velocity Track 


The second track that is created follows the outline of a line centered on the moon’s fixed axis 
latitude/longitude point (0 °N. 0 °E.). The length of the straight track is 30 km on the lunar surface, which 
would be the farthest distance that the rover would be away from the lunar module before it returns. This 
is approximately the distances that the Apollo rovers traversed in the later missions. The duration of the 
roving mission is set at 6 hr, 3 hr for each traverse of the track. The starting point of the track is placed at 
the point (0 °N. -0.496 °E.) (-0.496 °E. corresponds to 15 km on the lunar surface from 0 °E.). 

This straight track is designed to have a variable velocity throughout the simulation. Since the track 
has a length of 30 km and it takes 3 hr to traverse the track in one direction, the average velocity due east 
is 2.7778 m/s, but the track has a maximum velocity of 4.3771 m/s due east, due to the sinusoidal 
accelerations that are placed throughout the track. 

It is important to understand that this variable velocity track starts with zero velocity, and ends with a 
zero velocity when it reaches its destination, as a normal rover application would. While the rover is 
traversing the track in each direction, it is constantly accelerating or decelerating. Therefore, the rover 
does not have a constant magnitude velocity at any point of the simulation. In fact, the magnitude of the 
velocity profile follows the shape of the first half of a sinusoid waveform. 

Figure 2 provides an illustration of the straight track, which is created in Satellite Tool Kit (STK). 
Note that the simulation does not include any of the terrain information that is portrayed in the figure. 
Figure 2 attempts to illustrate the variable magnitude velocity nature of this track by capturing images of 
the track at 0, 0.75, 1.5, 2.25, 3, 3.75, 4.5, 5.25, and 6 hr after the starting point of the rover at the lunar 
module. Note that in the image bottom right region, representing the rover 6 hr after departure from the 
lunar module it returns to the starting point. 





Aag 




Aogi 






Figure 2. — Straight track velocity illustration. 
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3. Methodology 

This analysis uses the two rover tracks discussed earlier with IMU sensor data to propagate errors 
throughout the roving mission to determine the total error at the end of the rover mission. The error 
propagation analysis is updated every simulated second throughout the 6 hr roving application simulation 
(ref. 2). There are 20 distinct noise histories for each track that are propagated using Matlab. At the end of 
simulation for each track, plots are created showing the error growth in the three local topocentric 
coordinates. 

The basis for the error propagation simulation was derived in (ref. 2). Inputs to the simulation model 
include acceleration, case-to-inertial direction cosine matrix, and angular velocity. Outputs of the 
simulation are the same quantities corrupted by IMU errors. IMU errors can include the following: 

Gyroscope bias drift 
Gyroscope scale factor 
Gyroscope non- orthogonality 
Gyroscope random walk 
Accelerometer bias repeatability 
Accelerometer scale factor 
Accelerometer non-orthogonality 


Scale factor errors are the ratio of the input-to-output mapping for the gyroscope and accelerometer. 
Non-orthogonality errors arise because the actual gyroscope/accelerometer axes cannot be installed 
perfectly orthogonal to each other. In the gyroscope measurement Y g in the simulation, the IMU has an 
additive error 8 consisting of the following: 


Y g - 1 wf + 8 / wf 


( 1 ) 


8 ^wf -dc +S g J wf + r/wf 


+ Sv 


( 2 ) 


where 

1 wf true angular velocity 

dc gyroscope bias drift 

S g gyroscope scale factor matrix 

T g gyroscope non- orthogonality matrix 

s„, gyroscope Gaussian white noise random walk 

In the accelerometer Y a measurement in the simulation, the IMU has an additive error Sac consisting 
of the following: 

Y a -a c + 8a c (3) 

8a c =b a +Sa a c+r a ac+Ba a c ( 4 ) 
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where 


a c true acceleration 

b a accelerometer bias 

S a accelerometer scale factor matrix 

r„ accelerometer non-orthogonality matrix 

B a gyroscope-to-accelerometer misalignment matrix 

Now, if some of the parameters in equations (2) and (4) are rewritten as follows, equations (5) to (19), 
then equations (2) and (4) become equations (20) and (21). It is useful to have the parameters in these 
forms when computing the covariance of the IMU system. The equations for the covariance of the 
gyroscope and accelerometer errors are given in equations (22) and (23). 


d( ! w£)=diag( ! w£) 

s g ~ 

Sgw£ =D{ I wC)s g 
D{ac)-diag(a c ) 




ax,S ay fiaz 


S a a c = D(ac)s a 



IvuC 


w, 


0 

0 


0 

0 


Yg - \Fgyz, ^gyzXgzx, \ 

r g r wC = F^ w c)y g 


F i a c) = 


0 0 0 


-a x 0 0 


Cl X Cly 


Y a ~ }ayz, 'azyyazx, . 


(5) 

( 6 ) 

(7) 

( 8 ) 
(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 
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r a a c =F(a c )y a (16) 

0 a z -a v 

a*= — a z 0 a x (17) 

a v - a x 0 

K=[ B aX , B a y, B az\ (18) 

B a a c = -a*fi a (19) 

8 J wf = d c +D( I w^)s g +F( I wf)yg+s w (20) 

ha c = b a +D(ac)s a +F{ac}la-aZ$a ( 21 ) 

Pbwc^ 1 +P dc + D( 1 wf]p sg D T ( I wf)+F( 1 w^)+P yg F T ( I w^') (22) 

P5ac = Pba +D{da c )P sa D T {da c )+F{da c )P ya F T {5a c )+a^P tia af (23) 


where 

Pswc gyroscope error 

P^ c angular error in case coordinates 

Pdc gyroscope bias drift 

P sg gyroscope scale factor 

P lg gyroscope of the non-orthogonality 

P 5ac accelerometer error 

P ha accelerometer bias 

P sa accelerometer scale factor 

Py a accelerometer non-orthogonality 

Pp a gyroscope-to-accelerometer misalignment 

4. IMU Parameters 

The analysis performed in this document utilized the properties of an IMU model provided by 
James Russell Carpenter of NASA Goddard Space Flight Center (ref. 2). This IMU model contains 
gyroscopes and accelerometers. The best-case properties of the gyroscopes can be found in table 1, while 
the best-case properties of the accelerometers can be found in table 2. Additional parameters that were 
assumed for the simulations can be found in table 3. 
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TABLE 1.— GYROSCOPE BEST-CASE PROPERTIES 


Standard deviation of bias drift 

0.01 °/hr 

Standard deviation of scale factor 

33 ppm 

Power spectral density of random walk 

0.01 °/V(hr) 


TABLE 2.— ACCELEROMETER BEST-CASE PROPERTIES 


Standard deviation of bias repeatability 

0.06 milli-g 

Standard deviation of scale factor 

100 ppm 


TABLE 3.— ADDITIONAL PROPERTIES ASSUMED FOR ANALYSIS 


Standard deviation of gyroscope non-orthogonality 

0 ppm 

Standard deviation of accelerometer non-orthogonality 

0 ppm 

Standard deviation of initial angular error 

0 arc sec 


5. Results 

Circular Constant Velocity Track 

The resulting error graph for the circular constant velocity track is shown in figure 3. There are three 
subplots in the figure, representing the topocentric x, y, and z dimensions as the rover traverses the track. 
Since the rover is known to be on the surface, the topocentric errors could also be converted to latitude 
and longitude errors. This would eliminate the error in altitude if one could utilize a detailed elevation 
map of the lunar surface. Reducing this one error could reduce the overall error that is experienced in a 
simulation run. 

The unit for the error in any of the topocentric dimensions is kilometers. The maximum error 
experienced in any case for the circular constant magnitude of the velocity track in any of the dimensions 
is less than 4.5 km in the topocentric x dimension, 5.5 km in the topocentric y dimension, and 5 km in the 
topocentric z dimension. 


Straight Variable Velocity Track 

The resulting error graph for the straight variable velocity track is shown in figure 4. Again, there are 
three subplots in the figure, representing the topocentric x, y, and z dimensions as the rover traverses the 
track. Since the rover is known to be on the surface, the topocentric errors could also be converted to 
latitude and longitude errors. This would eliminate the error in altitude if one could utilize a detailed 
elevation map of the lunar surface. Reducing this one error could reduce the overall error that is 
experienced in a simulation run. The maximum error experienced in any case for the straight variable 
velocity track in any of the dimensions is less than 5.5 km in the topocentric x dimension, 4.5 km in the 
topocentric y dimension, and 4 km in the topocentric z dimension. 


NASA/TM— 2006-214354 


7 





Figure 3. — Circular constant velocity results. 



Figure 4. — Straight variable velocity results. 


6. Conclusions 

This analysis illustrates the error propagation that can occur on a lunar surface roving application 
when using an IMU that contains accelerometer and gyroscope biases and scale factor errors. This study 
is a first attempt to try and understand how these errors can affect performance of an IMU navigation 
system by creating two roving tracks with different traverse profiles and different magnitude velocity 
profiles. The study showed that it is possible to have errors up to 5 km in a single topocentric dimension. 
It is believed that the similar error propagation profiles are due to the nature of the rover traversing along 
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the spherical moon. When a similar analysis is done for lunar descent profiles to the surface, they have 
much different error propagation profiles (ref. 3). Using a combination of radiometric/IMU navigation 
system could reduce the error that is propagated, by closing the open-loop nature of the IMU system 
alone. Radiometric data would act as the feedback to make the radiometric/IMU navigation system a 
closed-loop system. 

There are a few types of actions that a follow on study should contain to advance the work performed 
in this study, such as: 

1. Create additional velocity profiles for these two track traverse profiles. This could include an 
acceleration/deceleration from a constant velocity. 

2. Simulate additional IMU performance characteristics to determine which error sources are the 
most dominant. 

3. Simulate these rover tracks with radiometric data from multiple lunar constellations to compare to 
IMU data alone. 

4. Simulate these rover tracks with a combination of IMU data and radiometric data to compare to 
IMU data alone and radiometric data alone. 
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